Head-on collisions of black holes: the particle limit 
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We compute gravitational radiation waveforms, spectra and energies for a point particle of 
mass mo falling from rest at radius ro into a Schwarzschild hole of mass M. This radiation is 
found to lowest order in (mo /M) with the use of a Laplace transform. In contrast with numerical 
relativity results for head-on collisions of equal-mass holes, the radiated energy is found not to be a 
monotonically increasing function of initial separation; there is a local radiated-energy maximum at 
ro ~ 4.5M. The present results, along with results for infall from infinity, provide a complete catalog 
of waveforms and spectra for particle infall. We give a representative sample from that catalog and 
an interesting observation: Unlike the simple spectra for other head-on collisions (either of particle 
and hole, or of equal mass holes) the spectra for 00 > ro >~ 5Af show a series of evenly spaced 
bumps. A simple explanation is given for this. Lastly, our energy vs. ro results are compared with 
• approximation methods used elsewhere, for small and for large initial separation. 
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I. INTRODUCTION AND OVERVIEW 



Among the earliest calculations of astrophysical sources of gravitational radiation is the "DRPP" calculation jlj (|] 
of radiation emitted when a particle, starting from rest at infinity, falls into a nonspinning black hole. In this work, 
the mass of the particle is treated as a perturbation of the Schwarzschild spacetime and the Einstein equations, to first 

a single linear wave equation. This work was 



s ! order in this perturbation, are organized into the Zerilli equation 

later elaborated by Ruffini and Ferrari J^] who considered the head-on plunge of a particle which starts with nonzero 
inward velocity at spatial infinity. 

There is now a renewed interest in such calculations. The advent of laser interferometric gravity wave detectors Q 
has directed attention at black hole collisions as one of the most plausible, and surely the most fascinating, sources 
of detectable waves. But the generation of strong waves is a process without simplifying symmetries and one which 
involves strong field interactions. The challenge of finding energies and waveforms generated in these processes helped 
q-( to spur the recent effort in numerical relativity and, in particular, in the Binary Black Hole supercomputing Grand 
Challenge @. Because of the difficulty of direct supercomputer solutions, approximate methods, as checks and guides, 
are very valuable. The "particle limit" in which the mass of one of the holes is very small, is in this category, and 
• • . models of holes in binary orbits are a fundamental tool. 

We are interested here in using a particle model to help in the understanding of numerical relativity results for 
collisions of holes. To date those computations have been limited to head on, axisymmetric collisions [||||. The codes 
have also been limited in the time for which they can evolve solutions. Infall from very large distances cannot at 
present be evolved; computations start from finite separation of the holes, so direct comparison with the DRPP results 
cannot be made. We present here an extension of the DRPP calculation to the case of infall from finite radius. 

The extension is not at all straightforward. If the particle starts at infinity, the initial value data can be taken as 
zero, so that the spacetime is unperturbed Schwarzschild until the gravitational influence of the particle is felt. If the 
particle starts from finite radius, the specification of compatible initial value data cannot be avoided. Understanding 
the influence of the initial data is, in fact, a major motivation for the present work. Supercomputer evolution will be 
limited, for the foreseeable future, to evolving the interaction of holes for short times. The starting point for these 
evolutions will be initial data from some approximation method (e.g., a post-Newtonian expansion). The way to make 
the bridge from the approximation to the strong field realm of supercomputing requires that we understand what the 
crucial features are of the initial value data. It seems to us that the particle limit provides a very useful probe of this 
issue. 

The remainder of the paper is organized as follows. In Sec. II we start by presenting the basic mathematical 
approach. Like the DRPP calculation, our approach uses frequency components. A Laplace transform converts the 
linearized Einstein equation for each £-pole, from a partial differential equation in Schwarzschild coordinates r, t to a 
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set of ordinary differential equations in r. The Laplace transform treats the initial data mathematically as if it were 
a source, and helps with insight into the relation of the initial value data and the stress-energy source. Two points 
in particular should be noted in Sec. II. The first concerns the linearized wave equation. The equation presented by 
Zerilli was formulated in the Regge- Wheeler gauge M|, a specific first order coordinate choice. The Zerilli equation 
was presented through a Fourier transform and has no direct equivalent as a partial differential equation. This is not 
a fundamental obstacle, but it does complicate the relationship of the Zerilli wave function and the initial data. To 
avoid these problems we use the Moncrief formulation of the problem Moncrief's wave function is constructed 
explicitly from initial data and is equivalent to the Zerilli wave function except for source-like terms. The second point 
to notice has to do with conventions. A variety of conventions has been used for normalizing the Zerilli (or Moncrief) 
functions and their transforms. In Sec. II, we present a definition of all our conventions and their relationship to those 
used elsewhere. 

Aspects of the numerical approach, a Green function integral, are discussed in Sec. III. A major issue is the need 
for highly accurate evaluation of the Green function integrals in order to preclude numerical errors from giving rise to 
an erroneous initial burst of radiation. Numerical results are presented in Sec. IV for energy, spectra, and waveforms. 
Plots are given for infall starting from rest at radii ranging from ro/2M = 1.1 to tq/2M = 15. These are combined 
with results for infall from infinity to give a complete picture of the dependence of the radiation on the initial physical 
conditions. The numerical results are discussed in Sec. V and are compared with predictions based on other results 
and approximations, especially those of the close-limit approximation. 

Throughout the rest of the paper we use c = G = 1 units, the signature h H — h, and other conventions of the 

Misner, Thorne and Wheeler textbook pj. 



II. MATHEMATICAL APPROACH 



A. Variables and conventions 



For the straight-line plunge of a particle into a nonspinning hole the perturbations of the Schwarzschild spacetime 
are all even-parity. We take the particle to be moving along the z axis, so that the perturbations arc not dependent on 
the azimuthal angle <fi, To describe the perturbations we use the notation of Regge and Wheeler |l0| , here specialized 
to the axisymmetric £-pole case (m = 0): 

ds 2 = ds 2 Q + (1 - 2M/r)(H Y m )dt 2 + (1 - 2M/r)- 1 {H 2 Y m )dr 2 
+r 2 (KY m + Gd 2 Y m /d8 2 )d8 2 + r 2 (sin 2 9KY m + G sin 6 cos 6 dY m / 86) d(f> 2 

+2H 1 Y m dtdr + 2h Q (dYm/d9)dtdd + 2hi(dY eo / 89)drd9 (2.1) 

where ds 2 , is the unperturbed line element for a Schwarzschild spacetime of mass M, where Hq, Hi, H 2 , ho, h\, K and 
G are functions of r, t, and where Yio(9) are the m = spherical harmonics. We then follow the prescription given by 
Moncrief (J] in his eqs. (5.8) - (5.12) and (5.27) to arrive at a wave function. Rather than choose the normalization 
of Moncrief's "Q" we choose a normalization closely related to Zerilli's [|| wave function. In the Regge- Wheeler |lC| ] 
gauge (G = ha = hi = 0) this "Moncrief- Zerilli" function is 

*< r '*> = ATT 

where we have used Zerilli's notation 



K 



r - 2M 
\r + 3M 



{H 2 - rdK/dr} 



(2.2) 



A= (£ + 2)(i- l)/2 . (2.3) 

We now give the relationship of this variable, which we shall use throughout the paper, to the following wavefunctions: 
(i) Q appearing in Moncrief's [Q eq. (5.27) (ii) -0pcrt and Vnum appearing in Ref. |l2[ . For 1 — 2 the wave function 
appearing as tjj in eq. (11-31) of Cunningham, Price and Moncrief [ pH] , and tjj defined in eq. (13) of Price and Pullin 
agree with V'pcrt- The relations are: 

^ = 2Q/[£(£ + l)] 

= 2[(*-2)!/(* + 2)!]Vwt ( 2 - 4 ) 
= v/2[(*-2)!/(* + 2)!]7/; num ■ 
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Zerilli || derives his equation only in the frequency domain, so no direct comparison can be made with functions of 
r,t. It is, however, straightforward to use the same steps as Zerilli in the Rcgge- Wheeler gauge, but confined to the 
r, t domain. It turns out that to decouple the even-parity variables in the perturbed Einstein equations, one needs to 
take an extra time derivative. One can then define the variable ipzeifr, t) by 



V'Zcr = 



r 2 K - (r - 2M)H l 
Xr + 3M 



(2.5) 



where the dot indicates a derivative with respect to t. The Fourier transform of this variable, divided by — iuj is the 
function defined in Zerilli's papers. For vacuum perturbations ^zer = V': so the Zerilli wave function, in vacuum, can 
be chosen to agree with our tp. 

By computing the stress-energy pseudotensor in a radiation gauge, we find that radiated power per solid angle is 



dPower 
dQ 



1 

64^ 



dtp 
du 



d d 2 

oote w Y eo - w Y io 



(where u is retarded time), so that the gravitational wave power, integrated over all angles, is 



Power 



1 (1 + 2)! 



dtp 



du 



64tt (I - 2)! 

Our wave function tp satisfies a second-order wave equation with a source term, 

d 2 tp d 2 tp 



dt 2 ^ dr* 2 



V e (r)tP = S e (r,t) 



Here r* = r + 2M\n(r/2M - 1) is the Regge- Wheeler (U| "tortoise" coordinate, V e is the Zerilli potential 

2M\ 2A 2 (A + l)r 3 + 6A 2 Mr 2 + 18AAf 2 r + 18M 3 



V t (r) 



1 



r 3 (Ar + 3M) 5 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



and Se (r, t) is the source term. 

We are interested in the case in which the source is a point particle, of mass mo, following a spacetime trajectory 
given by t = T(t), and r — R{t), where t is proper time along the particle world line, 



= (ma/U )U>*U v 6(r - R)S 2 (n)/r 2 



(2.10) 



Here = dx^/dr is the particle 4-velocity, and the 0-component is U° — Sq/(1 — 2M/r), where Eq = —Uq, the 
"energy-at-infinity per unit particle mass," is a constant of motion that takes the value \/l — 2M/tq for infall from 
rest at rg. The two dimensional delta function S 2 (fl) gives the angular location of the particle trajectory 



s 2 (n) = Ytm(0, <P)Yl m {e, & = Y m (eW(2e + i)/iir 



(2.11) 



with the last expression applying for infall along the positive z axis. 

We now n eed to write the the Einstein (G^) and Ricci (R^) tensors, in a tensor spherical harmonic decomposition 
like that in (2.1). The components we shall need are: 



Goo = GUY m Roo = J2 RUYm Ree = J2( R22KYm + R22Gd 2 Y m /d6 2 ) 



(2.12) 



(Here, and below, we omit I indices where the meaning of the symbols - such as G44 - is clear from context.) By 
rearranging the expressions for the perturbed Einstein and Ricci components, we find that the source term in fl2.S|) , 
for an axisymmetric problem, is 



S e (r,t) 



2(1 - 2M/r) 
r(A + l)(Ar + 3M) 

-(Xr + M)R22K- 



r 2 I 1 



2M 
r 



d 



R22K 



r 4 A 
Xr + 3M 



dr 

GU + r 3 RU 



(2.13) 
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The Einstein equations can now be used to replace the perturbed Ricci and Einstein components with components of 
the particle stress energy 



G44 = kU°(1 - 2M/r) 2 S(r - R)/r 2 
RAA = k [(1 - 2M/r) 2 U° - ^ (1 - 2M/r)] S(r - R)/r 2 
R22K = k[8{t- R)/2U°] , 



where 



k = 8Trm 0V / {2£+l)/ATT . 



(2.14) 



(2.15) 



When these are used in (2.13) the result is 

2(1 - 2M/r)n 



St(r,t) 



r{\ + l)(Ar + 3M) 



-r 2 (\-2M/r)—5'{r-R) 



2U° 



r(A + 1) - M 3MU°r(l - 2M/rf 



2U° 



rX + 3M 



8(r - R) 



(2.16) 



When the source term is included, the relation between our wave function ip, and the Zcrilli wave function of (2.5) 
becomes 



1p = IpZer 



KU°(r-2M) dR r/ 

— — 8(r - R) . 



(A + l)(Ar + 3M) dt 



(2.17) 



B. Initial data 

For comparison with established results of numerical relativity we want initial data representing an initially sta- 
tionary spacetime, so we take V'Oi the initial time derivative of ip, to be zero. To determine our ipQ, the initial value 
of tp, we choose an initial 3-geometry with the conformally flat form that is used in numerical relativity []l5| , [l6|| , 

ds 2 = $(f, Of [df 2 + f 2 (d0 2 + sin 2 9d(f> 2 )] , (2.18) 

where $ satisfies a flat-space Laplace equation V 2 $ = 0. We choose $ to represent a throat of mass mi on the z axis 
at Z = Z\ and a throat of mass mi at z — z% . The simplest such solution has the form 

+ , m2 1 . ,2.1!)) 




r 2 sin 2 9 + (r cos 9 — z\) 2 \j r 2 sin 2 9 + (r cos# — z-i) 2 



We now identify M = mi+m^, and choose Z\ — —z ijmi /mi) so that the dipole moment vanishes for f > z 2 - We treat 
the mass ratio m 2 /M as small and keep terms in ( [2.19 ) only to first order in this ratio. We ignore the nonradiative 
£ — 0, 1 perturbations and write $ as 

<£ = ! + -§ + 5 E m)Pe(™s6) (2.20) 

£=2,3,.. . 

where 

Mr) = \ i / -fl + i V- >Z2 • (2-21) 
w [ {r/z 2 ) + if r < z 2 K ' 

We next change radial variables, from isotropic-like f to a Schwarzschild-like coordinate r, with the transformation 

f = L/r + Vr - 2~m\ 2 /4 . (2.22) 
When this is used in ( [2.18 ), and terms higher order in m 2 are omitted, the result is 
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ds 2 = 1 



2m,2/r 
1 + M/2f 



Fe{r)Pe{cost 



e=2,3, 



dr 2 



1 - 2M/r 



+ r 2 (d9 2 + sin 2 6d<p 2 ) 



(2.23) 



From this we can infer that the initial value perturbtations in the Regge- Wheeler notation of (2T) are G = hi = 
and 



1 + M/2f" 



21+1 



(2.24) 



We must now complete the identification of this conformally flat solution with the point particle solution. Clearly 
the location f = Zi of the perturbative throat must be set to 



fo = [^+^r -2M^ /4 



(2.25) 



The mass m-2 of the throat must not be set to the mass mo of the particle. The particle mass mo can be viewed 
as the "bare" mass of the particle, the mass that is measured very close to the location of the part icle. On the 
other hand, mi is a formal parameter of the initial value solution; since the total ADM mass for (2.19) is mi + 7712, 
gravitational binding energy is included in TO2. We give here two arguments for the correct relationship. The fi rst 
is based on the definition, given by Brill and Lindquist of "bare" mass for a solution with the form ( 2.19| ), a 
mass exclusive of binding energy. To find an expression for this bare mass they look at the extension of the geometry 
through the perturbative throat into an asymptotically flat universe; the bare mass of the perturbative throat is the 
mass measured at infinity in that universe. From eq. (13) of Ref. |L5j the bare mass of the perturbative throat (to 
first order in mi/M) is 



m^ lc = m 2 (l + M/2f ) 
We identify the bare mass with the particle mass mo and conclude 



m 2 



mo 1 + 



M 
2fn 




(2.26) 



(2.27) 



An independent way of finding the relationship is to look at the perturbed hamiltonian constraint. For a conformally 
flat (Hq = H2 = K) initial 3-geometry [equivalent to eq. (C7a) of Ref. Jl7|] this gives us 



1 - 



2M\ d 2 K 



dr 2 



I 2M \ (2 3M \ 19K d 2M\£(£+1) 
r J r dr 



K = -GAA 



-87TTO0 



21+1 

47T 



U° 1 



r 

2MV 



T S(r - ro) 



(2.28) 



By integrating across the discontinuity, and by using j2~2l|),(p^), and ( gjg ), we find 

AK, r = dK/dr\ r=r+ - dKjdr\ r=r - 



-87TTO9 



V(21+1)/4tt 

-1/2 3/2 



(2.29) 



-87rmo 



T$y/l-2M/r 



which gives us the same relationship as in (2.27). This second derivation of the mass relationship is instructive. If 
we choose a conformally flat initial geometry then all the information about the perturbations of that geometry is 
contained in the single function K{r). But ( 2.28 ) completely fixes K(r), so the initial geometry has no freedom. 
(Alternatively: To choose different initial data we would have to introduce nonconformally-flat terms.) This is to be 
contrasted with the case of nonperturbative throats, where a variety of choices is possible for conformally flat initial 
data. 
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C. Laplace transforms 



We now define the Laplace transform \P of ip to be 

*(r,w) : 



^(r, t)dt 



(2.30) 



We take ip to vanish for t < 0, which means that \&(r, w) must be analytic in the upper half of the complex lu plane. 
At large r, for outgoing waves, tp(r, t) is a function only of t — r* so ^ takes the form 



(2.31) 



We call A the amplitude of the outgoing radiation. Since the outgoing radiation -0 (r - , t) is pure real, the amplitude 
satisfies the crossing relation A{—uj) — A*(uo). 

The waveform for outgoing radiation, as a function of retarded time u = t — r* can be found by the inverse transform 



+oo 



A(w)e~ 



"doj = Re 



A(w)e~ 



'dcu 



(2.32) 



For a particle falling in from infinity the waveform extends to u — > — oo, and there is no initial data to deal with. In 
this case A(uS) is to be interpreted as a Fourier, rather than Laplace, transform. If the particle has £o > 1, then the 
waveform does not vanish at u — > — oo, and the transform ex ists o nly for lu in the lower half plane, and will have a 
pole at u) = 0. For the inverse transform (the first integral in ( 2.32 ) the contour should be interpreted as going below 
the real lu axis, so the integral is the equivalent of the Cauchy principal value plus half the contribution of the residue 
of the pole at lu — 0, 



tp(u) = Re 



lim — 

e->0 7T 



A(w)e~ 



+ - lim ,H(w)} 

Z LJ—*0 



(2.33) 



From Parseval's theorem and (2.7), we have that the radiated energy at Scri+ is 

i 1 1 — ' ' i ' ( ■ / 1 1 1 ■ \ 

Energy : 



1 (1 + 2)! 
64tt (£-2)\ J_ 
1 (£ + 2)' r°° 



d S\ dt 



128tt 2 {i-2)\ 
1 (£ + 2)! 



du / 
uj 2 \A(lu)\ 2 dcu 



64tt 2 (£-2)! J Q 
The energy spectrum is therefore given by 

effinergy 1 (£ + 2)\ 



uj 2 \A{uj)\ 2 du 



64tt 2 (£ - 2) 



io 2 \A(u,)f 



(2.34) 



(2.35) 



We now multiply ( |2.8| ) by e tuJt , integrate from t = to oo, and integrate by parts in the integral involving d 2 ip/dt 2 . 

ial equation for 5*: 

2 + [lu 2 - V e (r)] * = -Mr) + twVto(r) + 5(r, w) , (2.36) 



The result is an ordinary differential equation for 5* 

d 2 V 



dr* 2 

where V'o(^) is the initial value of ^(i, r), and V'o('") is the initial value of ^>(t, r), and the source term S is defined as 

e^St(r,t)dt . (2.37) 



S(r,w) 

For the particle source we substitute ( [2.16 ) into ( 2.37 ), to get 

2(1 - 2M/r)K 



S{r,uj) 



r(A + l)(Ar + 3M) 



-r^(l - 2M/r) 



1 d 



|i?| \2U°R / 



1 



r(A + 1) - M 3MC7°r(l - 2M/r) 2 



2C/° 



rA + 3M 



(2.38) 



T(r) 
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The subscript 'T(r)" indicates that functions of time, such as R = dR/dt, U°, e luJt are to be evaluated at the value of 
t = T(r) where T = T(r) with r = R(t). For infall from rest at radius Tq, this gives i as a function of r through 



1/2 



r 
''() 



4M S / r \3/2 

+ Hm?J £ ° arctan 



+2 arctanh 



2M 2M 
r r 



(2.39) 



The result in ( 2.38 ) for S(r,uj), however, is not valid for r — ro. When r = r p the delta functions in ( 2.16 ), as 
functions of t, have vanishing arguments at t = 0, the endpoint of t integration in (2.37), so the integration is not well 
defined. The way to deal with this will be explained in the next section. 



III. COMPUTATIONAL IMPLEMENTATION 



A. Green function formal solution 



We start by rewriting ( 2.36| ) in the form 

a 2 * 

dr* 2 



+ [tj 2 -V e (r)] # = S tot (r,w) 



(3.1) 



where Stot is the complete right hand side of ( 2.36| ), including both the stress energy term and the initial value terms, 

Stot(w, r) = -ip (r) + iojifr (r) + S(r, u) . (3.2) 

This equation is to be solved for the boundary conditions of ingoing waves at the horizon, and outgoing waves at 
spatial infinity: ^ — * e~ lur for r* — > — oo, and $ — > e lur for r* — > +oo. The Green function solution is found in the 



usual way. We define yL{r*,uj) and y^(r*,uj) as the homogeneous solutions of (3.1) with asymptotic forms 



r — >oo 



yL(?-*,w) 

y R (r*,uj) — > e - 

We define the Wronskian of the homogeneous solutions, an r*-independent constant, to be 

W{u)) = vl-^Vr - VR-rrVL ■ 
dr* dr* 

With the above definitions, the Green function solution is written 
1 



W{lo) 

In the limit of large r* this gives us 



VR(r*,uj) Stot(r, uj)y L (r*, uj)dr* + y L (r* , w) / S to t(r,uj)y R (r* ,u)dr 



A(u>) = 



1 



S tot (r,u)y L (r*,u))dr* 



(3.3) 



(3.4) 



(3.5) 



(3.6) 
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B. Evaluation of the Green function integral 



The first two terms in ( |3.2[ ) are straightfor ward to integrate in (3.6), but the stress energy source term S(r,uj) 
cann ot be evaluated at tq. [The expression in ( 2.38| ) formally diver ges as (r — ro)~ 3 / 2 and hence cannot be used in 
( |3 .6| ) .] To make sense of this we start by writing the source in ( 2.1 6| ) as 



S(r, t) = F{r, t)S'(r - R[t])) + G{r, t)S{r - R[t])) 



(3.7) 



where F and G contain no delta functions. The troublesome part of the Green function integral can then be written 
as 



S(r,uj)y L (r* ,uj) dr* = X X + X 2 , 

where X\ is the integral involving G and involves F. The first of these is 

2 l = / +oo e <wt dt/^ o » L (r*,w)dr*G(r,t)*(r - R) 
= J °° e*"dt J 2 °° M VL(r*,uj)G(r, t)S(r - R)j^ 
= jo°° e^ t y L (r*(t),cj)G(R(t), t) dt/(l - 2M/R(t)) 
= - rL e<WT(r) ^( r * ,oj)G(r,T(r))dr* /R . 



(3.8) 



(3.9) 



In the final integral, the factor 1/R~ (r — ro) x ^ 2 diverges but is integrable. A similar set of transformations is now 
applied to X 2 , 



X 2 = 



7/f 



y L (r*,u)F(r,t)5'(r-R)dr* 



e lult dt / y L (r*,u)F{r,t)8'(r- R) 



dr 



2AI 



1 - 2M/r 



e~*dt 

Je>T{r) 







y L (r*(t),uj) — 



F 



F 







dr \l- 2M/r J ' (1 - 2M/r) 2 dr* 

F 



\ r J dr V 1 - 2M/r / 



y L (r*,uj) + 



d 



(3.10) 



r=R{t) 



1 - 2M/r J dr 

To evaluate these explicitly we need the fact that for free fall from rest at ro 



— {y L (r*,u)) 



dr* 
~R 



R 



2M\ 2M/r-2M/r 
' / 1 - 2M/r 



(3.11) 



We can now use the explicit expressions for F(r,t), G(r,t) from ( p. 16 ) to write 



Zi +z 2 



e iuT(r) dr * ^ _ 2M/r) 



dr* ^ L 



A + 1 y/2M/r - 2M/r (Ar + 3M) 

f\ M (2A-3 + 12M/r )Ml 

+^{a+i-- + ^ — Xr + J j 



From (3.6), (p.2j), and the definitions of X\ + X2, we have that 

1 



A{u) 



W{uj) 



Xi+X 2 + iuj ip (r)y L (r*,uj)dr* 



4>o(r)y L (r*,uj)dr* 



(3.12) 



(3.13) 



Since we are considering an initially stationary problem we have tpo = 0. The initial data for ip comes from putting 
( pT24[ ) in {$L~ 
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"00 



2m 2 ^4tt/(2£+1) 



A + 1 1 + M/2r Xr + 3M 



(A + l)r + M - rV 1 - 2M/r 



M/2r 



r e /f 



1 + M/2f J \ r e /f e +1 



(3.14) 



-n l 



2M 



-(£ + l)f£/r 



where the upper expressions apply in the case f > fo, and the lower for f < fo. 

It sh ould be noted that ipo approaches a nonzero constant, and yi(r* ,oj) — » e~ %u}r , as r — > 2M, so the integral over 
ipo, in (3.13) is improper at large negative r* . We must recall that we are reall y com puting A in the upper half of the 
complex uj plane. To deal with this computationally we recast the integral in ( 3.13| ) into the form 



M r )VL{r*,uj)dr* -> iuj- 1 ^ {2M)e- lur '^ 



4>o(r)yL(r* ,uj) dr* 



(3.15) 



where ip$(2M) is the limit of ipo at r = 2M. The value of r* tart must be chosen such that there is negligible variation 
of ip between r* tart and the horizon. 



C. Numerical method 



The hrst step in the solution is to determine the Wronskian in (3.4). We denote the form of y^ at large r by 



II L 



(3.16) 



We hnd /3(uj) by solving ( |3.l|) with the righthand side set to zero, and with the starting condition = e~ iwr imposed 
at a large negative value of r*. A fourth-order Runge-Kutta routine is used to integrate y^ out to large values of 
r* where it is matched to approximate forms of the asymptotic solution. In practice, good accuracy was difficult to 
achieve with the asymptotic form in (3.16) and asymptotic solutions one order higher in 1/ujr were used. From (3, 
the Wronskian follows immediately: 



W{lo) = 2iwP(w) . 

With (3 in hand, with ipo — 0, and with the substitution in ( p. 15 ), the problem consists of computing 



A( W ) 



2iu/3(uj) 



li+l 2 - A(2M)e^ ,ur ^" + iw 



ipo(r)yL(r* ,w) dr* 



(3.17) 



(3.18) 



where X\ + X2 is the integral given in ( pHt ) and V>o is given in fl3.14| ). A numerical sol ution for A(u>) is found by using 
a fourth-order Runge-Kutta routine to solve for y^ and dy^/dr*, and the integral in (3.18) is done by Simpson's rule. 
From the solution for A the energy spectrum is computed with (2.35), and the waveform from ( 2.32| ). 

The numerical solution used a routine to find the "particle contribution," X\ +I2, and one for the "initial value 
contribution," the integral over tp in (3.18). For both routines, second-order convergence was found and Richardson 
extrapolation was used. The step size in the Runge-Kutta and integration routines were halved until the Richardson 
extrapolate agreed with that from the next larger grid within a preset error limit. The initial value contribution could 
be usually be found within an error of 0.2%, while the particle contributions required an error preset of 0.5%. For 
most values of £, uj, tq, these precision requirements were easily met. The exception was a relatively small number of 
points at which the real and imaginary parts of a contribution differed by more than an order of magnitude. In this 
case it was difficult to get high accuracy in the smaller part. 

An estimate of the error in our results is complicated by the fact that the physically important results are a 
superposition of the particle and initial value contributions, and significant cancellations occur in this superposition. 
These cancellations, in principle, mean that the error may be much larger than the small relative error in each 
contribution. To arrive at an estimate of the error in our determinations of the radiated energy we have recomputed 
the energy for four trials, in which ±0.5% was added to the particle contribution, and then to the initial value 
contribution. The results for energies and errors are given in Table I. We see that for many cases the computed energy 
is not highly sensitive to the cancellation; the estimated error of around 1% is just what we would expect in the square 
(energy) of a quantity (amplitude) with an error of 0.5%. For some cases, however, especially those with higher £, 







there is a significant magnification of error. We emphasize that the error estimates given in Table I, are extremely 
conservative. In arriving at them we have used the maximum 0.5% error applied to all values of u>, whereas this 
maximum error actually applies only to a small subset of the points. The smoothness of (most) waveforms reported 
in the next section, and the consistent variation of results with changing ro, is evidence that the actual errors are 
rather smaller than those reported in Table I. 

In addition to the infall from finite ro we also have computed spectra and waveforms for infall from infinity. We 
characterize these cases with the same parameter sq we use for infall from ro. Here it has the value of the Lorcntz 
7-factor 



so 



i/Vi 



(3.19) 



for a particle with velocity Wqo at infinity. For computation of infall from infinity the above computational scheme is 
modified only in the following ways: (i) The form of T(r) in (2.39) must be changed to 



T(r) 



so 



: r i 1 (l'.u) 

(2e§ - 3)e 



el 



2M 
r 



hi 



1+ M- 1 »(a?) 



+ ln 



{el- 1)3/2 

~(2eg - 1) (r/2M) + 1 + 2e (r/2M) y/e 2 Q -l + 2M/r 



(r/2M) - 1 



(3.20) 



(ii) The initial value contributions in (3.18) must be omitted, (iii) The limit of integration in Q3.12 ) must be changed 
from rg to oo. The computed energy for these Sq > 1 cases are given in Table II. Since there is only a particle 
contribution in these cases, there is no issue of cancellation affecting the errors. The errors in the energy are, in 
fact, primarily due to the cutoff in the solution at a finite radius (with an analytic addition to represent the source 
contribution to infinity). Error estimates were made by varying the cutoff radius, and were found to be around 1% 
for all eq > 1 cases. 



IV. WAVEFORMS AND SPECTRA 



Results for quadrupole waveforms and spectra naturally divide themselves into three ranges, small ro/2M (less 
than ~ 2), moderate tq/2M (from ~ 2 to ~ 5) and large tq/2M. Waveforms are given as functions of retarded time 
u = t — r* . For small ro, as shown in Fig. [lja, the shape of the waveform is that of simple quasinormal ringing. This 
shape is the same for all ro and the single example shown suffices for all small ro- Since the waveforms have the same 
shape, the energy spectra, shown in Fig. [j]b, c, also have the same shape, changing only in magnitude as ro increases. 

As ro increases further, the early negative excursion of the waveform begins to broaden, and the spectrum shifts 
slightly. Fig. [2] shows the r$/2M = 2 spectrum along with the spectra for the two limiting cases, the DRPP infall 
from rest at ro — > 00, and the close-limit spectrum ro — > 0. (The latter is normalized to have the same energy as the 
ro = 2 spectrum.) The tq/2M = 2 spectrum has an appearance that interpolates between the two limits, as might be 
expected. As ro increases further, however, changes in the spectrum develop that might not be expected. As shown 
in Fig. ||, the simple spectrum for small ro develops a secondary peak and the secondary peak grows with ro. As ro 
continues to increase, the initial shape of the waveform becomes a very broad depression extending from the moment 
infall begins to the beginning of quasinormal ringing. This is illustrated in the waveform in Fig. ^a, for ro/2M = 10. 
The start of infall, at t = 0, r$/2M = 10 corresponds to u/2M = —12.2, and it is at this value of u that the waveform 
begins to take on nonzero values. The initial part of the waveform, then, represents the gravitational bremsstrahhmg 
from the early nonrelativistic part of the particle motion. The small wiggles around u/2M w — 12 are a numerical 
artifact due to imperfect cancellation of contributions from the initial value and particle parts of the source. (To 
verify this we changed the initial value contribution by ±10% and found the change in the initial wiggles to be much 
greater than in other features of the waveform). For comparison, a DRPP waveform, for infall from infinity, is also 
shown in Fig. ^. (For infall from infinity, of course, the zero of time cannot be set to the beginning of infall. The time 
was arbitrarily shifted for the DRPP curve.) This waveform has a similar ringing pattern as the infall from finite 
radius, but lacks the initial waveform depression. 

The end of the initial waveform depression is, roughly, the beginning of of quasinormal ringing, as can be seen in 
Fig. [|a. The generation of quasinormal ringing is associated with the peak of the potential in (|1]). The time 
t/2M rj 54 at which the particle reaches the peak, at around r/2M w 1.5 is at retarded time u/2M 53. This is 
consistent with Fig. 0a, which shows ringing beginning somewhere around this value of u. 
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It is interesting to compare the waveforms for a particle falling from rest, to the waveform of a particle on a time- 
symmetric geodesic trajectory, a particle that long in the past was moving radially outward just outside a hole, that 
reaches a certain maximum radius rg at ti me t = 0, and that subsequently falls into the hole. The analysis of this 



case requires only a simple modification of ( 3.15 ); the initial value terms are omitted and (due to time symmetry of 
the source) the complex conjugate of I\ + I2 is added to I\+l2- The resulting waveform is shown in Fig. ||a, and has 
also been given in Ref. . That waveform shows two periods of quasinormal ringing, an early one excited when the 
particle goes outward through the region around the potential peak, and a later one due to motion inward through 
the peak. This later period of ringing, and in fact all features of the waveform generated after u/2M w —10 agree 
very closely with the waveform for the particle falling from rest. Because the earlier ringing has a higher amplitude, 
the energy spectrum for the symmetric trajectory, shown in Fig. ^b, is dominated by this early ringing, and has very 
large amplitude. In the figure it is seen to be much larger than the DRPP spectrum, which (see Fig. ^) approximates 
the spectrum for infall from Tq/2M = 7.5. 

Since the spectra for infall from both small r$ and from 00 have a single peak (see, e.g., Fig. |2|) it is interesting that 
for intermediate values of r$ the spectra are characterized by a row of more-or-less evenly spaced bumps. As r$ grows, 
the number of bumps increases, and the bumps decrease in spacing and in height. The origin of these modulations 
can be understood by considering a simple example: Suppose a waveform consists of nothing but a single period of 
quasinormal ringing with the transform A(u>). The transform of the same waveform shifted later in time by T s hitt 
would be A((jj)e lujTshi{t . A wave consisting of two periods of ringing, the original and the later one, would then have a 
transform A(u>)(l + e lujT ' hitt ), and hence an energy spectrum 

w 2 |AH| 2 4cos 2 (^T shift ) . (4.1) 

The combined spectrum would have the shape of the single waveform spectrum modulated on a frequency scale 
Slu = 27r/T s hift- If the two waveforms were not identical, we would expect modulation of the spectrum, but less than 
100% modulation. 

This explanation can be tested on time-symmetric motion, where there are two well separated periods of ringing 
in in Fig. ||a. The amplitudes of the ringing are different, of course, so we shouldn't expect 100% modulation of 
the spectrum and, indeed, the modulations in Fig. ||b, are not 100%. The time shift between the first and second 
ringing periods is on the order of 75(2M) (roughly the time it takes for the particle to rise up from Tq/2M = 1.5 
to r /2M = 7.5 and fall back to r /2M = 1.5). This suggests that the spacing of the spectral bumps should be 
5uj = 27r/[75(2M)] = 0.84/2M, which is in good agreement with what is seen in Fig. ||b. The application to a 
nonsymmetric waveform, like that in Fig. [|a, is less obvious. The spacing of bumps 5ui « 0.11/2A/ in Fig. ||c suggests 
that the time shift is on the order of 57(2 M). This, presumably, represents the time between the descent of the 
waveform at u/2M w and the start of ringing at u/2M « 50. 

The interpretation works as well as a predictor of the bump spacing for other values of Tq, and seems clearly 
to be qualitatively correct, and to explain the progression of the spectra. As ro approaches the horizon, the time 
between the initial moment and the onset of ringing goes to zero, so the spacing of bumps is infinite and there are 
no modulations of the spectrum. We get the single-humped close-limit spectrum. As ro grows, the time between the 
initial motion and the start of ringing increases, so the spacing of the bumps gets smaller, and hence more bumps 
appear. As ro is becoming larger, however, the early waveform is becoming less dramatic (the initial depression is 
decreasing in amplitude) so in interacting with the later ringing it is producing smaller modulation; the height of the 
bumps is decreasing. Finally, as ro — > 00, we can think of the spectrum approaching one that has an infinite number 
of infinitesimally spaced zero-height bumps, bumps that are invisible in the DRPP spectrum. 

Results for higher ^-poles are shown in Figs. ^|. The waveforms for I > 2 show a more complicated structure of 
the pre-ringing radiation, resulting in more complex spectra. The total energy radiated in different ^-poles is shown 
in Fig. ^. Results are given for several different values of ro, and for particles falling in from infinity, both with no 
initial velocity (the DRPP case) and with £0 > 1. The distribution of energy among the multipoles is dominated by 
the quadrupole for small ro. As tq/2M increases to around 1.5, the higher multipoles become more important, but 
with further increase of r$/2M to ~ 2 the trend is reversed, and the ratio of multipole energies take on values that 
remain constant for further increases in ro. For infall from infinity, the ultrarelativistic cases, with high eq, radiate 
more heavily in higher multipoles, as would be expected. To find the total energy, given in Fig. we computed 
energy radiated in I = 2,3,4 modes and assumed that the energies from each £ decreased as a geometric series. This 
allowed us to add an estimate of the contributions from I > 4. The addition w as typically around 2% of the energy. 
The energy is plotted as a function of proper distance L — JJ*° drj yl — 2M/r, rather than ro, to show more clearly 
the details for small separation. 

Aspects of infall from infinity are given in Fig. ^[ which shows waveforms and spectra. The waveforms are char- 
acterized by a large early amplitude that is an increasing function of £o> and the spectra show strong low frequency 
radiation, extending to u> = 0, due to this early phase of the radiation. The total energy radiated in the lowest three 
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multipoles is also shown. (Since the relative importance of higher multipoles increases with sq, extrapolation to total 
radiated energy is not immediate.) The energy results are extended to include infall from finite radius where Eq takes 
the value — 2M/r . 



V. DISCUSSION 



One of the interesting questions that can be clarified with the above results is the validity of the close approximation. 
This approximation, for two holes, assumes that the holes are initially close enough so that the structure of the initial 
data at small radius is inside an initial all-encompassing horizon. Only the large r features of the initial data therefore 
are relevant to the production of outgo ing r adiati on. W e can immediately apply this method to the particle problem 



by comparing our initial geometry, in ( 2.23 ) and ( 2.21 ), with the initial geometry in equation (4.25) of Rcf 



12fl. We 



see that the results of that reference can be applied to the particle problem by the replacement 

8M^(/i ) -» m 2 (z 2 /M) £ . (5.1) 



With (2.25) and (2.27) this can be rewritten as 

1 ( mn \ / rn \ « / , — , \ 2 ^+! 



By the methods of Refs. 0,0], Eg/2M, the radiated energy in units of 2M, for each multipole, is shown to be 



1.26 x lCT 2 /^, 3.10 x 10 _3 k|, 8.33 x 10~ 4 k1, ■ -^respectively for I = 2,3,4- • •. If we replace n e with (|5_2j) we find 
the close limit predictions for energy. In Refs. Jl2],[l4| the outgoing radiation is computed from an evolution of ipi, 
using a finite difference representation of (|2.8|) (with no source term). Here we have also computed waveforms and 



energies in the close limit directly by the transform methods of Sees. II and III, with only the following changes: The 



integrals X\ + 2 2 are omitted from (3.13), and only the f > z 2 form of Ti(f~) is used in (2.21). The results of the 
two methods are energy values that agree to better than 1% and waveforms that are almost indistinguishable. (The 
numerically evolved waveforms are rather smoother than the waveforms from the transform method. They lack the 
small amplitude wiggles that can be seen, e.g., in the waveform in Fig. [l]a at early times.) 

In Fig. [l^, we plot the £ = 2 close-limit energy prediction and compare it to the full computation for the particle 
infall. It is clear that the close limit method is acceptable out to ro » 2.2M — 2.3M, and fails by a large factor at 
r = 3M. This is in accord with the general picture that the close limit sh ould be a reasonably good approximation 
when the particle starts inside the peak, around r « 3M, of the potential ( |2.9| ). The waveforms in the close limit all 
have precisely the same shape; only the amplitudes differ. The computed waveforms are in excellent agreement with 
the small ro waveforms (see, e.g., the waveforms for r$/2M = 1.5 in Fig. |l]a). The close-limit predictions of waveform 
shape are even better than the energy predictions. Only for ro/2M larger than around 1.5 does the waveform start 
to change from the close-limit shape. 

The dependence of radiated energy on ro, shown in Fig. §1 is perhaps the most interesting result of our computations. 
One has an intuitive instinct that the radiated energy should decrease with decreasing ro. At the crudest level this 
intuition is based on the idea that infall from a larger radius results in a particle which "strikes" the black hole harder 
and excites more quasinormal ringing. Since a decrease of energy with smaller ro is expected, a natural first guess 
for approximating the decrease is to multiply the DRPP energy by a reduction factor (1 — k2M jr®), where k is some 
fitting parameter of order unity. This approximation, with k = 1.5 is shown in Fig. [ll] and is compared with the 
computed results, and with the DRPP limit. The choice k = 1.5 was made to give good agreement at r§j2M =15, 
and presumably at larger values of ro- (It is difficult to compute energies at much larger values of ro, due to the rapid 
modulation of the spectrum.) 

A more physical justification for decrease of radiation with decreasing ro can be constructed starting with the 
quadrupole formula. A faster moving particle from a larger ro implies a larger value of the time derivatives of 
the quadrupole moment. This argument has been used |19[ | as the basis of a simple quantitative model for the 
effect of varying ro. The energy for infall from infinity, in that model, is reduced by a reduction factor F ro based 
on the quadrupole formula. [See eq. (23) of Ref. Pj|.] (For large ro that reduction factor reduces to F ro = 
1 - (60/27) (2M/r ) + O(2M/r ) 2 .) In Fig. |lT], we show the result of that simple model. It is clear that the F ro 
factor captures the correct qualitative feature of a decrease of radiation with decrease of ro but implies too dramatic 
a decrease. 

As ro continues to decrease, a rather unexpected effect appears. Below ro « 7M the energy begins to increase with 
decreasing separation. At yet smaller radii (ro less than around 4.5M) the energy again decreases with decreasing 
ro, as the close limit dictates it must. Thus the relationship of radiated energy and ro has the expected nature in 
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the two regimes where simple arguments apply: large separations and small separations. The anomalous behavior 
in the range A .5M -7M underscores the fact that the generation of outgoing radiation is tied closely to the nature of 
the potential ( |2.9[ ), which peaks around rg = 3M, and cannot be understood in terms of close or far approximations. 
This anomaly, it should be noted, appears to have no equivalent feature in the case of the head-on collision of two 
equal mass holes Jl2| . Presumably this is because the replacement of the particle by a hole means that the infalling 
hole is not localized at a particular value of the potential of the other hole. As the mass ratio of the infalling holes 
becomes smaller and smaller there must come a point at which an anomalous bump develops in the dependence of 
radiation on initial separation. 

We have seen that the particle limit provides a relatively easy tool for understanding some aspects of the generation 
of radiation, and of the collisions of holes. We intend next to use this formalism to study what features on initial data 
are important to determining how much energy is radiated for particle infall, and presumably for black hole collisions. 
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FIG. 1. The £ = 2 waveforms and spectra for small r . In (a) the waveform is shown as a function of retarded time 
u = t — r* , with t = corresponding to the moment at which the infall begins. Included for comparison is the predictions 
of the close-limit approximation. For the ro/2M = 1.1 case, shown in (b), the close-limit energy is larger than the computed 
energy by 11%. For the ro/2M = 1.5 case, in (c), the close-limit energy is 2.9 times the computed energy. 

FIG. 2. The £ = 2 spectra for ro/2M = 2, and for the two limiting cases: DRPP and close-limit. The close-limit case is 
normalized to have the same energy as the ro/2M = 2 spectrum. For the three spectra, close-limit, ro/2M = 2, and DRPP, 
the maximum of the spectra occur respectively at u> max (2M) = 0.765, 0.675, and 0.625. 

FIG. 3. Spectra of £ = 2 energy for ro = 3.5 and 4. For comparison, the DRPP spectrum for infall from infinity is also 
shown. 

FIG. 4. The 1 = 2 waveforms and spectra for large ro. In (a) the waveform for ro/2M = 10 is shown as a funtion of retarded 
time u, and is contrasted with the DRPP waveform (for which the zero of retarded time has a different meaning). Energy 
spectra are shown in (b)-(d) for three large values of ro and are contrasted with the ro — > oo DRPP limit (dashed curves) . 

FIG. 5. The 1 = 2 waveform and spectrum for a time-symmetric trajectory. In (a) the waveform (dotted curve) is given 
for a particle that moves outward and reaches a maximum at ro/2M = 7.5 before falling inward, and is compared with the 
waveform for infall from ro/2M = 7.5. In (b) the spectrum of energy generated by the time-symmetric motion is compared 
with the DRPP spectrum for infall from infinity. The DRPP spectrum contains less than 5% as much energy as that for the 
time-symmetric motion. 

FIG. 6. Results for higher multipole moments. For ro/2M = 5, waveforms and spectra are shown for £ = 2,3,4. Dashed 
curves are spectra for infall from infinity. 

FIG. 7. Energy in different multipoles. Energy for £ = 2, 3, 4 is shown for several values of ro in the case of infall from finite 
radius, and for several different values of e in the case of infall from infinity. The energies for ro/2M = 1.01 are multiplied by 
10 to improve the plot. 

FIG. 8. Total energy radiated by a falling particle, as a function of the initial proper distance of the particle from the horizon. 
The points shown are at r /2M = 1.01, 1.1, 1.2, 1.3, 1.5, 1.8, 2, 2.25, 2.5, 2.75, 3, 3.5, 4, 5, 6, 7.5, 10, 15. The local maximum 
is at ro/2M w 2.25 and the minimum at ro/2M w 3.5. 

FIG. 9. Results for infall from infinity. Plots are given of waveforms (a) and spectra (b), (c), for a particle falling in from 
infinity with nonzero energy. The total radiation emitted in the first three multipoles is shown as a function of eo- This curve 
also extends to particles falling from finite radius, with eo = (1 — 2M/ro) 1 ^ 2 . The local maximum of the spectrum occurs for 
0.75, and the local minimum at about so ~ 0.84. 

FIG. 10. Quadrupole energy for an infalling particle as a function of the particle's initial proper distance from the horizon. 
The computed energy is compared with the close-limit approximation. 



FIG. 11. Total radiated energy as a function of ro/2M. Computed results are compared with a simple model, and with the 
best l/r fit. 
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TABLE I. Radiated energy for infall from ro 
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